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Abstract. We study an iso-spectral deformation of general ma- 
trix which is a natural generalization of the Toda lattice equation. 
We prove the integrability of the deformation, and give an explicit 
formula for the solution to the initial value problem. The for- 
mula is obtained by generalizing the orthogonalization procedure 
of Szego. Based on the root spaces for simple Lie algebras, we con- 
sider several reductions of the hierarchy. These include not only 
the integrable systems studied by Bogoyavlensky and Kostant, but 
also their generalizations which were not known to be integrable be- 
fore. The behaviors of the solutions are also studied. Generically, 
there are two types of solutions, having either sorting property or 
blowing up to infinity in finite time. 
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1. Introduction 

In this paper we consider an iso-spectral deformation of an arbitrary 
diagonalizable matrix L £ 371 ( N, M) . With the deformation parameter 
i ER, this is defined by 

±L-P,q, (1.D 

where P is the generating matrix of the deformation, and is given by 
a projection of L, 

P = U(L):=(L) >Q -(L) <0 . (1.2) 
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Here (L) >0 , <0 , denotes the strictly upper (lower) triangular part of L. 
In terms of the standard basis of 9Jt(N, R), i.e. 

Eij = {8ikfije)i< k) e<N > ( 1 - 3 ) 
the matrices L and P are expressed as 

L = a >j E >) » ( L4 ) 

l<i,j<N 

P = a iJ E ij ~ £ • ( L5 ) 

l<i<j<N l<j<i<N 

Using the commutation relations for E^, i.e. 

, E ke \ = E ie Sjk — Ej k Sa , (1.6) 
the equations for the components = aij(t) are written in the form, 

dcij 



dt 




(1.7) 

where I := max(i,j) and J := min(i,j). The equation ( | . 1|) is also 
defined as the compatibility of the following linear equations with iso- 
spectral property of L, 

L$ = $A , (1.8) 

It® = P$) (L9) 
where $ is the eigenmatrix, and A is the diagonal matrix of eigenvalues, 

A = diag(Xi, ■ ■ ■ ,X N ) . (1.10) 

The set of equations (|1.8|) and (|1.9p is also referred as the inverse scat- 
tering transform for the system fll.lf) . The main purpose of this paper 
is to show the complete integrability of (1.1) with ( p..4|) and ( |1 . 5|) by 
means of the method of inverse scattering transform. 

One of the most famous and important example of Q1.1Q is of course 
the non-periodic Toda lattice equation, where L is given by a symmetric 
tridiagonal matrix |L3|. The matrices L and P for this equation are 
commonly written as 

N N-l 

L T = J2 a i E v + J2 b i(Ei,i+i + E i+1>i ) , (1.11) 

i=l i=l 
N-l 

P T = E k(E hi+1 -E t+1 ^ . (1.12) 

i=l 
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The integrability of the Toda lattice equation has been shown by the 
inverse scattering method JL2] [|L| . In this paper, we call ( |1 . 1|) with 
( |1.2| ) the "generalized Toda equation" . 



Several generalizations of the Toda lattice equation have been con- 
sidered. In 0, Bogoyavlensky extended the equation based on the 
simple roots of semi-simple Lie algebra q, where L and P were given 
by 

r 

i=i aen 
Pb = E & «(ea-e_ Q ). (1.14) 

Here the elements hi, e a , e_ a are Cartan-Weyl bases in g with r = 
rank(g) and II = the set of the simple roots. All of these equations 
associated with semi-simple Lie algebras are shown to be completely 
integrable hamiltonian systems. In JlT| Kostant solved these by using 



the representation theory of semi-simple Lie algebras. In [1[, Bloch 
et al. showed that these systems can be also written as gradient flow 
equations on an adjoint orbit of compact Lie group. They then showed 
that the generic flow assumes the "sorting property" (or convexity). 
Here the sorting property means that L(t) — > A = diag (Ai, • • • , Xn) as 
t — > oo, with the eigenvalues being ordered by Ai > A2 > • • • > A^. 

There are also other types of extensions: One of them is to extend Lt 
in ( |1 . 1 1| ) to a full symmetric matrix. The corresponding system, which 
we call the "full symmetric Toda equation" , was shown by Deift et al. 
P to be also a complete integrable hamiltonian system. In M Kodama 
and McLaughlin solved the initial value problem of the corresponding 
inverse scattering problem (|1.8| ) and (|1.9|) , and obtained an explicite 
formula of the solution in a determinant form. They also showed the 
sorting property in the generic solution. As a slight generalization of 
the full symmetric Toda equation, Kodama and Ye |10J considered a 



system with symmetrizable matrix L, which is expressed as L^y — L$S 
with a full symmetric matrix L5 and a diagonal matrix S. A key 
feature of this system is that the eigenmatrix of L^y can be taken 
as an element of noncompact group of matrices, such as 0(p,q), and 
defines an indefinite metric in the eigenspace. The integrability was 
also shown by a similar way as in ||, and the general solution now 
assumes either sorting property or blowing up to infinity in finite time 
as a result of the indefinite metric. 
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In ||, Ercolani et al. proposed the equation ( |1.1| ) with matrices, 

N-1 

L H = J2 + E b i3 E H > ( L15 ) 

i=l l<j<i<N 

P H = -2(L H ) <0 = -2 E » ( L16 ) 

l<j<i<N 

which was called the "full Kostant-Toda equation". This is also an 



extention of the Toda equation (|1.11|) which can be written in the 
form, 

JV-l N N-1 

l T = E + E a ^« + E b l E i+hi ■ (1-17) 

i=l i=l i=l 



As we will show in this paper, the transformation from ( |1 . 1 1| ) to ( |1 . 1 7| ) 



is given by a rescaling of the eigenvectors of Lt- Several group theoret- 
ical structure of the extended system were found. However the ques- 
tion whether the system is completely integrable still remains open in 
a sense of explicit integration. 

It is immediate but important to observe that all of these exten- 
sions are special reductions of the generalized Toda equation ( p..l|) . In 
fact, we show that these reductions are obtained more systematically 
as crtain decompositions of the root spaces of semi-simple Lie algebras. 

In this paper we prove the complete integrability of any reductions 
of the generalized Toda equation ( |1 . 1| ) . The content of this paper is as 
follows: We start with a preliminary in Section 2 to give some back- 
ground information necessary for analysis of the system ( |1 . 1| ) and the 
inverse scattering scheme ( |1.8| ) and ( |1.9| ). 

In Section 3, we solve the initial value problem of ( |1.9[ ) for the general 
system ( [Lip by generalizing the method developed in || and PH| . A 
key in the method is to use the orthonormalization procedure of Szego, 
which is equivalent to the Gram-Schmidt orthogonalization method. 
This shows the complete integrability of the generalized Toda equation 
in the sense of inverse scattering transformation method. 

In Section 4, we present reductions of Ql.lQ according to the classi- 
fication of semi-simple Lie algebras. The matrix L here then contains 
"all" the root vectorts, and it gives a generalization of the system for- 
mulated by Bogoyavlensky || . A key ingredient here is to find a matrix 
representation of the algebra in a decomposition consisting of diago- 
nal, strictly upper and lower matrices (Lie's Theorem ||). Then the 
integrability of these systems associated with semi-simple Lie algebras 
is a direct consequence of the result in Section 3. 
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Section 5 provides other reductions which include the full Kostant- 
Toda equation and system with a matrix L having band structure in 
the elements. 

In Section 6, we discuss the behavior of the solutions. Generically, 
in addition to the sorting property, there are slutions blowing up to 
infinity in finite time, as in the case discussed in |T(J . 

Finally we illustrate the results obtained in this paper with explicit 
examples in Section 7. 



2. Preliminary 

Here we give some background information necessary for the inverse 
scattering method (|1.S| ) and ( |1.9|) . As we will see in the next section, 
a key idea for solving these equations is to use an orthogonality of the 
eigenfunctions of ( |1.8|) . This is simply to consider a dual system of 
(|1.8|) and (|1.9| ), which are written by 

= , (2.1) 

= _ (2.2) 

dt 

where the matrix \1/ is taken to be and of course 

= = (2.3) 

In terms of the eigenvectors, these matrices can be expressed as 

$ = [</>(Ai), ••• ,<K\ N )] = M^)]i< itj < N • (2-4) 

* = [V(Ai), •■• ,V(Aiv)] T = [ipj(K)]i< itj < N ■ (2.5) 

Note here that 0(Aj) and are the column and row eigenvectors, 

respectively. Then the equations (|2.3|) give 



N 

k=l 
N 

E0i(A*)^i(A fc ) = 5^ , (2.7) 

k=l 

which are called the "first and second orthogonality conditions" . With 
(|2~7j), one can define an inner product < •, • > for two functions / and 
g of A as 

TV 

</,<?>:=E/(A*)<?(A fe ), (2-8) 
fc=i 
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which we hereafter write as < fg >. From L = $A\^, the entries of L 
are then expressed by 

a i:j := (L) y =< \<\>ii\)j > . (2.9) 

This gives a key equation for the inverse problem where we compute L 
from the eigenmatrix $ (and \l/) with the eigenvalues Aj. So the eigen- 
matrix plays the role of the scattering data in the inverse scattering 
method. Then the method for solving the initial value problem of the 
hierarchy (|1.1| ) can be formulated as follows: First we solve the eigen- 
value (or scattering) problem (|1.S| ) at t — 0, and find the scattering 
data, $° := $(0). Then we solve the time evolution of the eigenmatrix 
from (|1.9| ), and with the solution we obtain L(t) thorough the 
equation (|2.9|) . 

3. Inverse scattering method 

In this section, we construct an explicit solution formula for the 
initial value problem of the generalized Toda equation ( [Lip by using 
the inverse scattering method. A key of this method is to generalize the 
orthogonalization procedure of Szego with respect to the inner product 
(|2~8|). This is essencially an extention of the method developed in [[|. 

Following U we first "gauge" transform $ and \I> by 

$ = G<f> , * = §G (3.1) 

where the matrix G is given by 

G = diag < 0i^i > -1 / 2 , • • • , < (j) N ip N >~ 1 / 2 

Note that the gauge transform ( p.l|) includes a freedom in the choice 
of and ip, that is, Q3.1Q is invariant under the scaling 0j , ^ — > 
fi(t)<fii, fi(t)ipi, with arbitrary functions of t. With (|3TT|) , the 

equations ( |1.8|) and (|1.9Q , as well as (|2.1|) and ( p.2|) , become 



(G-^LG) $ = $A , * (CLG" 1 ) = A^ , (3.2) 

s*=( G " IpG )*-(s kl8G )*- (33) 

|* = -4( G P G -)-*(| 1 „ gG ) . (3.4) 



Noting G 1 (L) <0 G = (G 1 LG) <0 etc, we write 

G- X PG = -2(G- 1 LG) <q + G- 1 LG - diag(L) 
GPG- 1 = 2(GLG- 1 ) -GLG- 1 + diag(L) , 
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from which we obtain the equations for the column vectors <j)(\,t) in 
<3? and the row vectors ip(\,t) in 

^ = -2 (G~ X LG) <q + A0 - (diag (L) + ^ log g) , (3.5) 



/f -2V> (GLG- 1 j >q + Xifj - if) [diag (L) + - log Gj . (3.6) 

We here observe that ( |3.5|) and ( |3.6|) can be split into the following sets 
of equations by fixing the gauge freedom in the determination of (f> and 
if). In the components, these are 



~ 2 2^ rj + x< Pi > I 3 - 7 ) 

i i 



!§l = -2 i f:*l^ft> + ^, (3.8) 

o* i= i < Mi > 

~ log < > = an . (3.9) 



It is easy to check that ( p.7|) or (fl.8|) implies ( |3.9|) . It is also immediate 
from (13.71) and (B.8I) that we have: 



Proposition 1. The solutions of ( \3. 7\ ) and ( ITS) can be written in the 
following forms of separation of variables, 

4>{\,t) = M{t)(fP{\)e xt , (3.10) 
f{X,t) = iP°(\)N(t)e xt , (3.11) 

where M(t) and N(t) are, respectively, lower and upper triangular ma- 
trices with diag[M(t)] = diag[N(t)} = I , the identity matrix. 

Note here that the initial data for </> and if) are chosen as those of <f> and 
ip, i.e. 0(A,O) = 0(A,O) := 0°(A) and if)(X, 0) = ^(A,0) := ip°(X). As 
a direct consequence of this proposition, and the orthogonality of the 
eigenvectors, (|2.7|) , i.e. < faifjj >= for i ^ j, we have: 



Corollary 1. (Orthogonality): For each i,j G {2, ■ ■ • ,N}, we have 
for all t e R, 

< 4>iip° e e 2Xt > = , for 1= 1,2, ••• ,i- 1 , (3.12) 

<0°^e 2A *> = , for k = 1, 2, • • • ,j — 1. (3.13) 

Now we obtain the formulae for the eigenvectors of L in terms of the 
initial data {$(X)}i<i<N and {i>j(X)}i<j<N- 
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Theorem 1. The solutions (f>i(\,t) and ipj(\,t) of ( \3. %) and ( \3. 8j ) are 

given by 



A-i(t) 



Dj-i(t) 



Cn ... Cij-i $_(\) 

cu ... Ci,i-i <Pi(X) 
c n ... cij 



Cj-l,l 

V>?(A) 



(3.14) 



(3.15) 



where Cij(t) =< (j^ip^e >, and Dk(t) is the determinant of the k x k 
matrix with entries Cjj(t), i.e., 



Cy(t) 



i<*,i<fc. 



(3.16) 



(Note here that Cj,(0) = 5jj and A(0) = l.j 



Proo/. From equation (gl2|) and ( gig) with ( plQ|) and ( gUP , we 
have 



£ M«(t) < 0^°e 2A * >= , /or 1 < < i — 1 , (3.17) 

£=1 

£ < 0^ fc °e 2A * > N kj (t) = , /or 1 < £ < j - 1 , (3.18) 



k=\ 



Solving ( P7T7D and (|3l8|) for M ie and iV fci with M a = % = 1, we 
obtain 

' 1)<+ ' n ^1 (*):=#4t> 1<^<*, (3-19) 



M u {t) 
N kj (t) 



A-i(t) 
(_l)fc+i 

A-iW 



A 
A 



A-i(t) 

(t) :- 



A-i(t) 



1 < A; < j ,(3.20) 



where A 



is the determinant of A after removing the £-th row 



and z-th column. From ( |3.10| ) and ( fj.ll ), we then have 



= 1 



fc=i 



A-i(t) fi 



5>?(A)A, fi (t) , (3.21) 



A-iW fc =! 



EWlW' (3-22) 



Noticing that A^j is the cofactor of the element c& of the matrix 
{c mn ) i<m,n<i an d Aj , fe is f or Cjfc of (c mn )i< min <j we confirm ( |3.21| ) and 
O) are just ( CT) and (|3~T5|) . □ 

We then note: 



Corollary 2. T/ie gauge factors < (fiiipi > can be expressed by 

Proof. From ( |3.21| ) and (|3.22|) , we have 

1 



(3.23) 



<MA)V*(A) 



n 2 

u i-i e,k=i 



£ A^A^(A)^°(A)e 



2Af 



Then taking the bracket <, > in ( |2.8| ) leads to 

1 



i,kC£k ■ 



n 2 

Using the relation S|=i A^c^ = Z^i^^ with A^j = A-i complete the 
proof. □ 

This yields the formulae for the normalized eigenfunctions 



<t>i{\t) 



-,\t 



A(t)A-i(t) 



Cn • 
<ki ■ 

Cll 



Cj — 1,1 

V>i°(A) 



c M -i 0?(A) 
Ci,i-i 0°(A) 

Cij 



€(A) 



(3.24) 



(3.25) 



With the formula (|3.24 ) and ( [3.25 ), we now have the solution (|2.9| ) of 
the inverse scattering problem ( |1.8|) and ( |1.9| ). 

The above derivation of the eigenvectors is the same as the orthog- 
onalization procedure of Szego ||16fl , which is equivalent to the Gram - 
Schmidt orthogonalization as observed in ||. 

Remark 1. The generalized Toda equation ( |1 . 1|) with ( |1.2|) possesses 
a hierarchy defined by 



^-L=[P n , L] , n = l,2, 



where P n is given by 



P n = U(L n ) ee (L n ) >0 - (L")<o 
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(3.26) 
(3.27) 



The commutativity of these flows can be shown by the "zero" curvature 
conditions of P n , i.e. 

BP BP 

W L ~W L + [P m ^ p n]=0, (3.28) 

which is a direct consequence of the choice of (|3.27Q [|i~0|1 . The solution 
for the hierarchy is then obtained by extending the argument At in the 
solution formula to f (A, t) := En=i PH- 

Remark 2. In ||14|| , [ p"5|j and ||, the authors considered the following 
generating matrix P for the equation ( [[ . 1|) with the general matrix 
LG<m(iV,M), 

P = (L) >0 - (L T ) <0 = (L) >0 - [(L) >0 ] T . (3.29) 

They then showed that this equation is integrable in the sense of a ma- 
trix factorization of QR type, and the solution has the sorting property. 
Our method developed in this section can be also applied to this prob- 
lem as follows: First we note that the product $*$ of the eigenmatrix 
$ and its adjoint $* := $ is invariant under this flow ( |1 . 1|) . Then we 
define a hermitian matrix S = {sij)i<ij<N a s the inverse of $*$, i.e. 

= / . (3.30) 

The matrix S is determined from the initial eigenmatrix $ , and S<&* 
gives the inverse of that is, we have S$>* for \1/ in our method. 
Note that if L is symmetric, S is an identity matrix / and $ G O(N). 
In general, we see from the Binet-Cauchy theorem that S is positive 
definite. The equation ( |3.30|) now gives the orthogonality relation, 



0i(Afc)sH0j(A^) = Sij , (3.31) 

l<k,£<N 

from which we define the following inner product as in 



</,</»:= E fih)skeg(Xe) = < 9 J > • (3.32) 

l<fe,£<iV 

This leads to a positive definite metric. Then following the procedure 
in this section with some modifications based on \1/ = 5$*, we obtain 
the same result for the eigenvectors (|3.24j ) except the quantities Cij 
which is now given by 

Cij =< $ e x \ 0°e A * >= Cji . (3.33) 

The solution L(t) is then given by L(t) = $AS*$*, i.e. 

a ij{t) <f>j > {t) ■ (3.34) 

Thus, we can show explicitly the integrabihty of the equation ( |1.1|) with 
the generator P given by (|3.29|) for arbitrary diagonal matrix L, and 
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as a result of the positivity in the metric, the solution has the sorting 
property. The detail will be discussed elsewhere. 



4. ISOSPECTRAL FLOWS ON SIMPLE LlE ALGEBRAS 

In this section, we consider the generalized Toda equations ( |1 . 1| ) 
associated with simple Lie algebras q, and show their integrability. 
The matrices L and P here are given by a generalization of Q1.13Q and 



i=l aGA+ (3eA~ 

p s = b ^ e » - H c P e p ■ ( 4 - 2 ) 

aeA+ /3eA- 

Here hi are the bases for the Cartan subalgebra with r = rank(o), 
A + and A~ are the sets of positive and negative roots with the corre- 
sponding root vectors e a and ep{= e_ a ). These vectors {hi,e a } form 
the Cartan- Weyl bases of the simple Lie algebra g which satisfy for 
i,j G {1, • ■ ■ ,r} and a, (3 G A := A+ U A - 

[hi, hj] = 0, [hi, e a ] = a(hi)e a , 

[e«, e p \ = N a pe a+(3 , if a + (3 E A , (4.3) 

fe„,e_J = h n , for a e A + . 



Using representations of the Cartan- Weyl bases, we now express ( 4.1 
and Q4.2D in matrix form for each simple Lie algebra. Then we prove 
that the equation (|1.1| ) with those L s and P s associated with the Lie 



algebra q is completely integrable by the inverse scattering method 
developed in Section 3. The key ingredient in the proof is to show that 
for each simple Lie algebra q there exists a "permutation" matrix O s 
such that the matrices L s and P g are similar to L and P in ( |1 . 1| ) with 
P defined by (|L"2|) , i.e. 

L = O L g Oj, (4.4) 

P = O g P g O T s = U(L) . (4.5) 

In another word, we look for a similarity transform such that the matrix 
representations of e a for a G A + and e@ for (3 G A - are transformed 
to strictly upper and lower triangular matrices, respectively. The exis- 
tence of such representations is due to Lie's theorem 0. Here we give 
its explicit formula for each semisimple Lie algebra. Then the result 
in Section 3 implies the integrability of the system ( |1 . 1| ) with L g and 
P s . Note here that the generalized Toda equation is invariant under 
the similarity transform. Here we consider all the classical simple Lie 
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algebras, A n , B n , C n and D n . The system associated with the excep- 
tional algebra can be treated as the same way. For convenient matrix 
representations of the Cartan-Weyl bases, we follow the notations in 
i and @. 

A n : Let Eij be the (n + 1) x {n + 1) matrix defined in (Of ). We 
then take an element of the Cartan subalgebra as h = J2i=i \Eu with 
E7=i \ = 0. Using (|rj) for E tj , we have 

[h, E tJ ] = (Aj - Aj)£Jy . (4.6) 

Thus gives a root vector corresponding to the root a(h) = Aj — A 3 -. 
The simple roots are defined as 

= A fc - Ajfc +1 , for k = 1, • • • , n. (4.7) 

Then the positive (negative) roots are given by Aj — Aj with i < j 
(i > j). This implies that the choice of the Pa h is the same as that in 
(|1 . 2| ) . Note also that adding some constant to the Cartan subalgebra, 
one can choose hi of the basis to be En. Namely, the generalized Toda 
equation ( [Lip with (|1.4|) and ( |1.5| ) can be considered as an iso-spectral 
flow on the simple Lie algebra A n . 

C m : The element of this algebra is given by a 2m x 2m matrix X 
satisfying X T J + JX = where J is defined by 

J=f_°? Im ). (4.8) 

Here m is the m x m 0-matrix, and I m is the m x m identity matrix. We 
then choose the following bases with the 2m x 2m matrix E^ defined 
in O, 

= — Ej +m> i +m , 1 < i, j < m , 
e - • = £?ij +m + E j)i+m , 1 < z < j < m , (4.9) 



-ij J -'%+m,3 i -^^+771,1) 

Writing h = YhLi \ e u as a general element of the Cartan subalgebra, 
we have 



[h,eU = (A<-Aj)eL i^j 



[Ke%] = (Ai + A^eJ, <<j, (4.10) 
[/i,e|.] = -(Ai + A^eJ, z<j . 
The simple roots are taken by 

oik{h) = Ajfc - Afc+i, for 1 < fc < m - 1 , (4.11) 
a m {h) = 2A m , 
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from which the sets of positive and negative root vectors £j and 
are given by 



J Crt 



J Crr 



{ 4'» 
{ e--» 



2 

iji ^kl 
3 

ij > C H 



1 < j < « < m, 1 < /c < ^ < m} 



(4.12) 
(4.13) 



Then the matrix Lc m can be represented by 

L c = 



A x A 2 
A 3 



(4.14) 



where A%, 



A<± are the m x m matrices satisfying the relations 



A\ = -A 4 , A 2 = A\, A 3 



Al 



The matrix Pc m is now given by 

PCm = 

We then obtain: 



n(A) 

-A 3 



A 2 

-n(A 4 ) 



(4.15) 
(4.16) 



Proposition 2. VFzi/i the permutation matrix Oc m , we have the gen- 
eralized Toda equation ( |i.i[ ) on C m with L-P pair given by 



L 
P 

where Oc m is given by 



O 



O c L c Oi 
O c P c O] 



n(L) 



r*. 



with the m x m matrix Q 

( 


Q 



0„ 

Qr. 



1 \ 

1 



(4.17) 
(4.18) 

(4.19) 



Q 



(4.20) 



: : 
VI ... 0/ 

Proof. From ( |4.14j ) and ( |4.16| ), it suffices to show 

- Q m Yi{A A )Q T m = n(Q m A 4 Q^) . (4.21) 

Recall that the multiplication of Q m from the left (right) implies an 
exchange of rows (columns). Then we see 

<3m(At)>0(<0)Qm = (Qm^4<5m)<0(>0) , (4.22) 

which implies the assertion. □ 

Note that the equation ( |1.1|) with L Cm and Pc m is just a reduction of 
the generalized Toda equation on A 2m -i- 
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Example 1: We take the simplest case C 2 . The matrices Lq 2 and Pc 2 
are represented as 



j c 2 



( cti 




&2 


64 \ 


Cl 


a 2 






c 2 


c 4 


—ax 


-Cl 


V c 4 


c 3 


-61 


-a 2 ) 



(4.23) 



and 



c> 



( foi 6 2 64 \ 

-ci 64 63 

-c 2 -c 4 ci 

V -c 4 -c 3 -61 y 



(4.24) 



Under the similarity transformation with Oc 2 defined in ( [4.19| ), Lc 2 
and Pc 2 becomes 



Oc 2 L C2 Oc 2 



and 



Oc 2 Pc 2 Oc 2 





( 0,1 


61 


64 


6 2 ^ 






Cl 




&3 


&4 






c 4 


c 3 


— a 2 


-61 






V c 2 


c 4 


-ci 


-ai y 




( ai 


&1 


64 




\ 




-Cl 


o 2 




64 






-c 4 


-C3 


— a 2 


-h 




\ 


-c 2 


-c 4 


Cl 


-ai / 



(4.25) 



(4.26) 



Note here that under the similarity transformation the root space is 
decomposed into the diagonal, upper and lower triangular parts of the 
matrix (Lie's theorem). 

D m : The matrix representation of this algebra is given by a 2m x 2m 
matrix X satisfying X T K + KX = 0, where K is defined by 



K 







m * rn 
I'm. Om 



(4.27) 



Then the bases can be chosen as 



e 2 - 



E{j — Ej +m ^ +m , I <i,j < m 



E 



i,j+m 



Ei 



i+m,j 



E jt i +m , 1 <i < j <m , 
E j+m ,i, 1 < i < j < m . 
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(4.28) 



With a general element h = J2iLi ^i e li °f the Cartan subalgebra, we 
have 

[h,e\ 3 ] = j*k, 

%e%] = (Ai + A^, i<j, (4.29) 

from which the simple roots may be taken as 

ock(h) = Afc — Afe+i, for 1 < k < m — 1 , (4.30) 
+ A 

The sets of positive and negative root vectors £j m and are then 
given by 

£+ m = {4,4 | 1 < i < j < m} , (4.31) 
^D m = {el p e%\l<j<i<m}, (4.32) 
The matrix Ljj m is then expressed as 



JDm \ A, A 



Ai A 2 

3 ^4 



(4.33) 



where the m x m matrices Ai, • • • , A4 satisfy 

Af = -A 4 , A 2 = -A^ A 3 = -A^. (4.34) 
The matrix Pn is 



n(Ax) A 2 

-a 3 -n(A 4 ) 



(4.35) 



It is then immediate to see from Proposition 2 that the permutation 
matrix Oo m is the same as in the case of C m . Namely we have: 

Proposition 3. With the permutation matrix Oo m = Oc m given in 
( U-1Q ), we have 

L = Dm L Dm O T Dm , (4.36) 

P = Dm P Dm O T Dm = U(L) . (4.37) 

B m : The element of this algebra satisfies the same relation as in D m , 
X T K + KX = 0, except now K is the (2m + 1) x (2m + 1) matrix 
defined by 

1 T T 

h = [ m I m I . ( 1.38) 

lui Of/ 
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where is the m-column vector with entries. This algebra is referred 
as the orthogonal algebra so(2m+l), while the algebra D m is as so(2m), 
and has the same bases as ( 4.28 ) with additional elements, 

el = E 0i - E i+mfi , 1 < i < m , 

e\ = E i0 - E QMm , 1 < i < m , (4.39) 

where we have labeled the indices of the matrix Eij as < i,j < 2m. 
With the expression h = YJILi \ e \i as in the case of D m , we have 

[h,4\ = -Xiet, [h ) e\] = \ l el 

The simple roots are then chosen as 

otk{h) = A fc -A fc +i, for l<k<m-l, 
a m {h) = X m . 

The sets of positive and negative root vectors are now 



(4.40) 
(4.41) 



{el 



I P 1 P 3 P 4 
I '-'iji '-■jit °fc 



1 < 2 < j < m, 1 < ^ < fn} 
1 < j < i < m, 1 < k < m} 



(4.42) 
(4.43) 



The matrix is then expressed as 

/ 

L>B m = 



°1 b 2 



b 2 At A 2 | , (4.44) 
V -&i A 3 A 4 

where 6i, 62 are the m-column vectors, and the mxm matrices A\, - ■ ■ , A4 
satisfy the same relations as (^4.34j) . The matrix Pg m is now given by 

-bl bl 



-62 U(A 1 ) 
bi -A 3 




(4.45) 



We then have: 



Proposition 4. With the (2m + 1) x (2m + 1) permutation matrix 
C>B m , we have 



L = 
P = 

where Os m is given by 
Ob. 



-Dm -Dm -Dm ' 



OnLnOi 



B Pb O 



U(L) . 




(4.46) 
(4.47) 



(4.48) 



Proof. Under the similarity transformation with 0# m , L and P in 
(|4.46|) and (|4.47|) are given by 

/ A l -b 2 A 2 Q m 
L= | b T x b T 2 Q m ^ | , (4.49) 
\ Qm^-s —Q m bi Q m A 4: Q rn 



and 



P 



U(A 1 ) 

lT 



A 2 Qr 



~b 2 

-b{ 

-Q m U{A 4 )Q 

III 



b 2 Q r 



(4.50) 



Then the equation fl4.21Q immediately leads to the result. □ 



Example 2: We take the simplest case B 2 , where Lb 2 and Pb 2 are 
represented by 



J B 2 



and 



b 2 



1 


c 3 


c 4 


-6 3 


-64 \ 


h 




6i 





b 2 




Cl 


a 2 


-b 2 





-c 3 





-c 2 


— CL\ 


-Cl 


\ -c 4 


c 2 





-h 


-«2 / 


/ o 


-c 3 


-c 4 


-h 


-64 \ 


h 





h 





&2 




Cl 





-h 





c 3 





c 2 





Cl 


V c 4 


-c 2 





-61 


/ 



(4.51) 



(4.52) 



Under the similarity transformation with Ob 2 defined in ( 4.48 ), L and 
P in flP9D and ( fcgCD become 



and 



(4.53) 



P 



/ ai 


61 


63 





-62 \ 


Cl 


a 2 


&4 


b 2 





c 3 


c 4 





-w 


-63 





c 2 


-c 4 


-a 2 


-61 


V -Pa 





-c 3 


-Cl 


-fli / 


/ 


61 


&3 





-62 \ 


-Cl 





&4 


b 2 





-c 3 




c 4 


-64 


-63 







c 2 c 4 





-h 


V c 2 





C3 


Cl 


/ 
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(4.54) 



5. Reductions on root spaces 



As we have explained in the introduction, several generalizations of 
the Toda equation may be obtained by taking reductions of the gener- 
alized Toda equation ( |1 . 1| ) with general matrix L. We then showed in 
the previous section that the equations on simple Lie algebras studied 
in are generalized by taking all the root vectors in the algebras. In 
this section, we consider reductions of these equations by restricting 
the set of roots in the sums in ( |4.1| ). 

Let S + and S~ be subsets of positive and negative roots of a simple 
Lie algebra q defined by, for Va £ S + and W/3 G S~ , 

S + := { a G A + | a -< a }, (5.1) 

S- := { (3 G A- | (3 y p }. (5.2) 

Here "-<" and 'V are the standard partial orderings between roots. 



We then consider the equation (LI) with the matrices L and P given 
by 

n 

L = Y a i h i+ Y b ^a + Y C /3 e /3> ( 5 - 3 ) 
i=l aeS+ /3e5" 

P = Y b a e * ~ Y C P e P> ( 5 - 4 ) 

where n = rank(o). We here claim: 

Proposition 5. The equation with L and P is a reduction of the 
generalized Toda equation on g. 

Proof. All we need to show is that the commutator [P, L] is in the 
span of the root vectors whose roots are in S + and S~. From ( p.3|) and 
(|5.4|) , [L,P] can be written as 

n 

[P,L] =-Y,a i [h i ,P) + 2 £ Y b ^[e a ,e p ] (5.5) 
»=1 aes+ /3e5- 

Using ( |4.3|) we first note that the terms [hi, P] does not produce any new 
root vectors. The second term, which then gives [e a ,ep] = N a pe a+ p, 
has a root a + (3 G S + U S~ (if a + (3 G A) , since a G A + and (3 G A~. 
This completes the proof. □ 

Example 3: The generalized Toda equation with band matrix L. 

This example can be obtained as the following reduction on A^-i- 

Consider the subsets of the roots S + and S~ given by 

S+ = { g A + | < j -i < M + < N- 1} , (5.6) 

S- = { g A - | < i- j < M~ < N- 1} , (5.7) 
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where M + and M are some positive integers. Then the corresponding 
matrix L which we denote L( M + M -\ becomes 



■>(M+,M-) 



( an 

a l+M~,l 



V o 



a l,l+M+ 







o>n,n-m- 



0-N-M+,N 



(5i 



O'NN 



J 



As a special case of this example, we now costruct the full Kostant- 
Toda equation having Lh-Ph pair given in ( 1.15 ) and ( p.. 16 ). Here we 
choose S + and S~ to be the sets of the simple roots (i.e. M + = 1) and 
of all the negative roots (i.e. M~ = N — 1), respectively. Thus the 
corresponding matrix is expressed as 



Then we claim: 



/ On h 

«21 0,22 

V a, m ■■■ 





b 2 











\ 



&N,N-1 a NN j 



(5.9) 



Proposition 6. The generalized Toda equation with L^jv-i) and 
P(i.n-i) '■= n(L( 1 jv_ 1 )) is equivalent to the full Kostant-Toda equation 
having Lh and Ph '■= —2(Lh)<o ^ n ( \l-13j) and \l.lb\ j. 



Proof. Let if be a matrix defined by 



H = diag 



1, h, &1&2, 



N-l 

n * 

i=l 



(5.10) 



Then the matrix Lh can be written as Lh = HL^i m-^H' 1 , which 
corresponds to the rescaling of the eigenvectors, i.e. = if $ with 
Lh§h = &hA. From the equation fll.Tp for L^^-i), we have 



^ = (a i+lji+1 - au)bi, for i = 1, 



N-l 



from which H satisfies 
d 



dt 



H 



diag(Ln) — aniiv )H . 
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(5.11) 



(5.12) 



Note here that diag(Ln) = diag(L^ 1)N _x^) = diag[au 
the derivative of Lh is calculated as 

' diag(L H ) , L H 



dt 



+ 



H 



(Inn]- Then 



(5.13) 



where we have used dLn N _^/dt = [P^n-i), £(i,jv-i)L an d HP(i,n-i)H 1 
U(Lh)- Noting the relation 

II(Ltf) + diag{L H ) = L H - 2{L H )<o , (5.14) 

we complete the proof. □ 

Thus the full Kostant-Toda equation can be solved through the gener- 
alized Toda equation with the £(i,jv-i)--P(i,jv-i) P arr as the reduction 
on Aff-i, that is, with the solution Ln jv-i), Lh = HLi\n-\\ 



'(1,7V-1), — 

6. Behaviors of the solutions 



Here we study the behavior of the solution of the generalized Toda 
equation obtained in Section 3 by following the approach in [1(1 . Many 



results obtained in |T0[ are valid for this more general situation. First 
we note: 

Lemma 1. The determinants Di for i = 1, 2, • • • , iV in are real 

functions. 

Proof. In the construction of the solutions $(t) and ^{t), the "gauge" 
G is fixed by (p|). In terms of A, @ is 



1 d Di 

log . 

2dt A-i 



(6.1) 



Note that D Q = 1, -Dj(O) = 1 and an are real functions. Then we see 
by induction that all Di are real functions. □ 

In (|6.1|) , suppose A (to) — for some finite t and some i. Then if L(t ) 
is a finite matrix, D^ito) must be also 0. By induction, Di(t Q ) = 0, 
but D(j(t) = 1, this forces an to be infinite, which is a contradiction. 
So we have: 

Lemma 2. Suppose Di(to) = for some to < oo and some i. Then 
L(t) blows up to infinity at to. 

We note that Di for i = 1, 2, ■ ■ • , N are the i-th leading principal 
minors of the product of matrices $ e \I/ e , where $ e and \l/ e are defined 
by 



/ e A ^0?(A 1 ) 



\ 
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e A ^*0?(Ajv 



and 



/ e Al V?(Ai) 



,Ait„/,0 



Vft(Ai) \ 



V 



e Ajv Vi°(A,v) ... e x ^%(X N ) J 
Then from the Binet-Cauchy theorem, we have: 



Lemma 3. The determinants Di with i 
as 



Di(t) =5> 2 £t=iV 

JiN 



1,2, 



N can be expressed 



(6.2) 



■01 (An] 
V>i°(A*) 



$(a*; 



where = (ji, ■ ■ • represents all possible combinations for 1 < 
ji < • • • < ji < N. In particular D (t) = 1, andD N (t) = exp(2^ili A^t) . 

This Lemma is very useful to study the asympototic behavior of Di for 
large t. 

We now obtain: 

Theorem 2. Let the eigenvalues of L be all real and ordered as \\ > 
X 2 > ■ ■■ > Ajv- Suppose that det(Q° k ) ^ and det(^° k ) ^ for k = 
1, . . . , N, where $° and are the k-th leading principal submatrices 
o/$° and respectively. Then ast^ oo, the eigenf unctions <j)i(\j,t) 
and ipj(Xi,t) satisfy 



(f)i(\j,t) -> S {j x 



det^det^U) 



Sij x 



deti^det^l^U) 

det^Jdetj^) 
det^^det^l^U) 



(6.3) 
(6.4) 



which implies the sorting property as t —>• oo ; that is, L(t) = &(t)Aty(t) 
A. 



Proof. Here we give a proof for (|6.3| ). The case for ipj(Xi,t) is ob- 
tained in the same way. Using Lemma 3, and from the ordering in the 
eigenvalues we see that the leading order term for Di is given by 



Di(t) -> e 2 ^=i Afc *det($°^°), as t -» oo . 
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(6.5) 



From ( 3.24 ) and ( p.5| ), the eigenfunctions behave as t — > oo 



<MA;t) 



The dominant term in the determinant gives 



#(A*J 



(6.6) 



di(t) •■• ci,i_i(t) 0?(A) 



Cii(t) ■■• Ci,i-i 0°(A) 



0?(A £l ) ... 0?(A, s _J 



0?(A) 



0?(A) 

2(6.7} 
#(A) 



Noting that the de- 



/ 1 2 • • • A; 

where P& is the permutation \ n B 

y ti 12 ■ ■ ■ tit j 

terminant in ( |6.7| ) is zero for A = Xj, with j — 1, . . . , % — 1, we complete 
the proof. □ 



This theorem implies that if all the eigenvalues of L are real, then 
generic solutions have the "sorting property" in the asymptotic sense. 
It should be however noted that Di(t) might be zero for some "finite" 
times, where the solution blows up (Lemma 2). Next theorem provides 
sufficient conditions for the solutions to blow up to infinity in finite 
time. 



Theorem 3. Suppose some eigenvalues of L are not real, det& n ^ 
and deffin ^ 0, for n = 1, • • • , N. Then L(t) blows up to infinity in 
finite time. 

Proof. We order the eigenvalues of L by their real parts. We still 
assume all the eigenvalues to be distinct. Since L is a real matrix, the 
complex eigenvalues appear as pairs. For a convenience, we also assume 
that there is at most one pair having the same real part. Suppose 
Afc + iflk and A^ — if3k are the first pair of complex eigenvalues. Then 
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= 2 Ej=i x i t + 2i ht 



from (|6.2| ), the leading order term in D k is 

• • 0i(A fe + i/3 fc ) 

• • <P° k {h + Wk) 
.. d>l{\ k -if3 k ) 

.. 4>l{\ k -i(3 h ) 
Since D k is real by Lemma 1, one can write the above as 



0?(Ai 



■01 (Ai 
€(Ai 

<(Ai; 
^(Ai; 



^(A fc + i&) 
V?(A fe -z/3 fc ) 



A cos(2/5 fc t) + 5 sin(2/? fc t) 



where A and B are two real constants. The above is an oscillating 
function about zero. Thus by Lemma 2, we conclude that L(t) blows 
up to infinity in finite time. □ 

7. Example 

In this section, we demonstrate the results obtained in this paper by 
taking an explicit form of the matrix L. The main purpose here is to 
solve the generalized Toda equation ( |1.1| ) for this explicit matrix, and 
discuss the behavior of the solution. 

Let us consider a 2 x 2 matrix L{t) = (<kj)i<i,j<2- The generalized 
Toda equation then gives 



d l 


a n 


Ol2 ' 




It \ 


k 021 


«22 J 





(7.1; 



The initial data of L(t) is assumed to be 

L(0) = 



1 

a b 



(7.2) 



where a and b are arbitrary constatnts. The eigenvalues of L(0), Ai 
and A 2 , are 

' (7.3) 



A 



1,2 



-{b±VV 



4a 



Then the initial eigenmatrices $° and \1/ are expressed by 



1 1 

Ai A2 



A2 — Ai 
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A 2 -1 
-Ax 1 



(7.4) 
(7.5) 



In order to compute the solutions $(t) and *&(t) from ( |3 . 24| ) and (|3.25|) , 
we need the quantities c^- =< (j) ip°e 2Xt >. From (|7.4|) and (|?T5|), they 
are 



Cll(t) = 


1 


A 2 — X\ 


ci 2 (t) = 


1 


A 2 — Ai 


C 2 l(t) = 


A1A2 


A2 — Ai 


c 22 (t) = 


1 


A2 — Ai 



A 2 e 2Al * - A ie 2A2 ' 

2A 2 t 



(7.6) 



e 2Al * - e 



A ie 2Al * + A 2 e 



2X 2 t 



from which the determinants -Dj(t) in (|3.16|) become 



D 1 {t) = c 11 (t), D 2 {t) 



cu(t) c 12 {t) 



C2\{t) C 22 (t) 

We now have the solutions (Theorem 1), 



3 2(Ai+A 2 )t 



$(*) 



e Al * 



D A 2 < 



Di(t) 



Aie Aa * A 2 e Al * 



(X 2 -Xi)JD 1 (t) 



A 2 e Al * 
-A ie A2i 



A 2 t 



-e 

D Aii 



(7.7) 

(7.8) 
(7.9) 



The solution L(t) of the generalized Toda equation is then obtained 
from O, Oij(t) =< A0iVj > (*), 



L(t) 



1 



A 2 e 2Al * - Aie 



2A 2 i 



AiA 2 (e 2Al * - e 2A2t ) (A 2 - A!) e ( Al+A2 )' 
-AiA 2 (A 2 - X 1 )e^+^ t A 2 e 2Ali - A 2 e 2Aai 



(7.10) 



Now let us discuss the solution behavior for t > 0. First we assume 
both eigenvalues Ai and A2 to be real. With the choice of the eigen- 
values in (|7.3|), we have Ai > A2. Then if A1A2 < 0, then the function 
Di(t) doesnot vanish for all t. This implies the sorting property (The- 
orem 2). For the case of Ai > A2 > 0, the D\ vanishes and we have the 
blowing up in the solution at the time t = tg > 0, 



f 1 A i 



(7.11) 



This formula also implies that for > Ai > A2 we have the sorting 
result for t > 0. Note here that the blowing up occurs at one time 
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t = ts (|7.11|) , and then the solution L(t) will be sorted as t — > oo, with 
the asymptotic forms of the eigenmatrices, i.e. ( |6.3| ) and (|6 . | ) . 

(7.12) 
(7.13) 

For the case of the complex eigenvalue Ai = A 2 := a + i/3, Di(t) is 
expressed as 

D^t) = e 2at sec 6 cos (2{3t + 6) (7.14) 

with tan# = a//3. This indicates the blowing up (Theorem 3). 

In the case of degenerate eigenvalues Ai = A 2 (i.e. b 2 + 4a = 0), we 
take the limit A 2 — > Ai := Ao in ( [7.1 Op , and obtain 

/ -2Agt 1 \ 

= r^A7 • (7 - 15) 

1 " 2A °* V " A o 2A (1 - A„*) J 

which showes the "sorting property" as t — > oo, i.e. L(t) — * Ao/ 2 . It 
should be noted however that L(0) with the degenerate eigenvalues is 
not similar to the "diagonal" matrix Ao/ 2 . 

We summarize the above results in Figure, where we classify the 
behavior of the solution in terms of the parameters a and b in ( |7.2| ). 

Put Figure 

In the figure, the shaded region corresponds to the blowing up solutions 
for t > where the eigenvalues are either complex or real with Ai > 
A 2 > 0. The other region including the positive b- axis and the lower 
boundary of the curve b 2 + 4a = gives the sorting property. 
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